On the Performance of 1 -level LDPC Lattices 



Mohammad-Reza Sadeghi 
Faculty of Mathematics and Computer Science, 
Amirkabir University of Technology 
Email: msadeghi@aut.ac.ir 



Amin Sakzad 
Department of ECSE, 
Monash University, Victoria, Australia 
E-mail: amin.sakzad@monash.edu 



o 

cn 



C/2 



> 
in 
o 

O 

m 



X 



Abstract — The low-density parity-check (LDPC) lattices per- 
form very well in high dimensions under generalized min-sum 
iterative decoding algorithm. In this work we focus on 1-level 
LDPC lattices. We show that these lattices are the same as lattices 
constructed based on Construction A and low-density lattice-code 
(LDLC) lattices. In spite of having slightly lower coding gain, 1- 
level regular LDPC lattices have remarkable performances. The 
lower complexity nature of the decoding algorithm for these type 
of lattices allows us to run it for higher dimensions easily. Our 
simulation results show that a 1-level LDPC lattice of size 10000 
can work as close as 1.1 dB at normalized error probability 
(NEP) of 10"'' .This can also be reported as 0.6 dB at symbol 
error rate (SER) of 10"* with sum-product algorithm. 

Index Terms — LDPC lattice, generalized min-sum algorithm, 
PEG algorithm, lattice decoder. 

I. Introduction 

P OLTYREV ifTTl suggests and investigates coding without 
restriction for infinite arrays such as lattices on the 
AWGN channel. That is a communication without power con- 
straints. In this communication scheme, instead of the coding 
rate and capacity two new concepts of normalized logarithmic 
density (NLD) and generalized capacity Coo are defined. For- 
ney et al. |9 1 showed the existence of sphere-bound-achieving 
and capacity-achieving lattices via Construction D lattices 
theoretically. He also established the concept of volume-to- 
noise (VNR) ratio as a parameter for measuring the efficiency 
of lattices. Therefore, generalized capacity for lattices means: 
there exists a lattice of high enough dimension n that enables 
transmission with arbitrary small error probability whenever 
VNR approaches 1. In addition, it has been shown in ifTTI 
that, this error probability is bounded away from zero when 
VNR < 1. A capacity-achieving lattice can raise to a capacity- 
achieving lattice code by selecting a proper shaping region Q, 

El. 

The search for sphere-bound-achieving and capacity- 
achieving lattices and lattice codes has begun with [12]. LDPC 
lattices are those that have sparse parity check matrices. These 
lattices were introduced first by Sadeghi et al. [12|. In this 
class of lattices, a set of nested LDPC codes are used to 
generate lattices with sparse parity check matrices. Another 
class of LDPC lattices, so-called LDLC lattices introduced and 
investigated in IS), IITtI . Turbo lattices employed Construction 
D along with Turbo codes to achieve capacity gains flSl, fT6^. 

An extended version of Progressive Edge Graph (PEG) 
algorithm, called extended PEG (E-PEG) algorithm llT2l can 
be employed to construct the Tanner graph of regular LDPC 



lattices. Since presence of short cycles can hurt the per- 
formance of lattices llT4l . such construction provides us the 
Tanner graphs for lattices with high girth. This will have an 
impact on the excellent performance of LDPC lattices. 

In the present work, we introduce and investigate 1-level 
regular LDPC lattices fT2l. We show that these lattices are 
equivalent to lattices constructed based on Construction A. We 
also provide a relationship between these lattices and the well- 
known LDLC lattices. Finally, we provide experimental results 
to reveal the effectiveness of 1-level regular LDPC lattices in 
terms of fundamental coding gain and error probability. 

II. Preliminaries on Lattices 

In order to make this work self-contained, general notations 
and basic definitions are given next. 

A. Lattices 

A discrete, additive subgroup A of the m-dimensional real 
space is called a lattice. We denote the length of the 
shortest nonzero vector of A by d^i^{K). Indeed, d^i^{K) 
refers to the minimum distance between different points of the 
lattice. Every lattice A has a basis B — {bi, . . . , b„} C R™ 
where every x G A can be represented as an integer linear 
combination of vectors in B. Let span(A) be 
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then a Voronoi cell v{-x.) can be defined as the set of those 
points of span(A) that are at least as close to the point x as 
to any other lattice point. The set of all vectors in span (A) 
where their inner product with all the vectors of A is in Z 
forms another lattice, denoted by A*, called the dual lattice 
of A. Define the coding gain of A as: 
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(det(A))^' 

where the det(A) is the volume of the lattice. The parameter 
det(A)~ refers to the normalized volume of an n-dimensional 
lattice A [9|. This volume may be regarded as the volume 
of A per two dimensions. Suppose that the points of a lattice 
A are sent over an AWGN channel without restrictions with 
noise variance . The parameter 

det(A)t 



VNR = 



(2) 



is the volume-to-noise ratio (VNR) of lattice A. When n is 
large enough, the parameter VNR can be interpreted as the 



ratio of the normalized volume of the lattice to the normalized 
volume of a noise sphere of squared radius ncr^. This is also 
defined as generalized signal-to-noise ratio (SNR) in 1 12] and 
in g]. 

Let the vector x in A be transmitted on the unconstrained 
AWGN channel, then the received vector r can be written 
as r = X + e where e = (ei,...,e„) is in the Euclidean 
space and its components are independently and identically 
distributed (i.i.d.) Gaussian random variables with zero mean 
and variance . Thus the probability of correct decoding is 
given by 

1 [ -iit|p 
Pc = / e^^dt, (3) 

where ||x|| is the Euclidean norm of x. Due to the geometric 
uniformity of lattices and without loss of generality we sup- 
pose that the vector is sent and the vector r is received. Then 
the components of r are i.i.d. Gaussian distributed random 
variables with zero mean and variance a^. The probability of 
error under a maximum likelihood decoder for A is the same 
as the probability that a white Gaussian n-tuple r with noise 
variance fall outside the Voronoi region i/(0) — v, i.e. 
Pe — 1 — Pc- The normalized error probability (NEP), P*, 
of an n-dimensional lattice A is the error probabiUty per two 
dimensions, i.e. P* = f Pe 13, ifTSl . 

B. Geometric structure of lattices 

We suppose that the n-dimensional lattice A' is a sublattice 
of an n-dimensional lattice A, and the lattice A' has a 
basis along the orthogonal subspaces S = {Wi}"^! with 
dim(yVi) = 1. For every point x in A, the label sequence 
is .g(x) = (51 (x), . . . ,.9„(x)) where .gi(x) = Pw^W + Aw, 
and (A) is the projection onto the vector space Wi and the 
cross section Ayy- is AnWi. The label code of A is then the set 
of all possible label sequences £ — g{A) = {g{x) : x e A}. 
It is a group code over the alphabet Q = Z^j^ x • • • x Zg^ , 
i.e. £ C Q. For every a.,c E Q, one can define ||2| the inner 
product (a, c) by: 

(a,c) :-^ + --- + ^ (modZ) (4) 

.91 9n 

where the operations like multiplications and divisions are 
done in M. For every subgroup £ ^ G, the dual of £, is 
defined by £* = {c G ^ | (c, a) ==0, Va e £}. Based on the 
above discussion, over the coordinate system S = {Wi}f=i, 
a lattice A can be decomposed as IS), lfT4l 

a = z"c(A) + /:p(A) (5) 

where /! is a group code over "^gi X 

X • • • X Zg^ and 
C(A) = diag(det(Awi), ■ • . ,det(Aw„)) 
P(A) = diag(det(Pwi (A)), . . . , det(Pw„ (A))) 

where diag(- • • ) denotes a diagonal matrix. The decomposi- 
tion (|5]l means that a vector v e K" is in A if and only if 
it can be written as v = zC(A) + cP(A) for some z e Z" 



and c e £. It is also shown ||2l that A* can be represented in 
terms of £*, P(A) and C(A) as follows: 

A* Z"P-1(A) + £*C-i(A). (6) 

In other words, let £ be the label code of A, then the dual of 
£, denoted by £*, is the label code of A*. It has to be noted 
that g*, the i-th alphabet size of £*, equals gi, the ith alphabet 
size of £, for 1 < i < n. 

The above representations for A and A* have a practical 
consequence that we can derive a generating set of lattice 
A* based on a generating set for £* and vice versa. More 
specifically, assume that {v^, . . . , v* } be a generating set for 
A*, then a generating set {c|, . . . , c*} can be derived easily. 
To obtain the j-th coordinate of c*, it is sufficient to divide 
the jth coordinate of v* by (j, j)-th entry of C^^(A) which 
is 1/ det{A\Yj) and then evaluate the result module gj. On the 
other hand, if {c*, . . . ,c*} be a generating set for £*, then 
£*C^^(A) along with the following vectors form a generating 
set of A*: 

( OTTTAJ). 0, 0, 0, 0) 

( 0, 0, 0, 0) 

( 0, 0, 0, 0, dct(Pv!.„(A)) )- 

To construct a Tanner graph for group code £ |2|, we use 
the following set of check equation constraints, which fully 
describes the codewords c — {ci, . . . , c„): 

y e z, l<k<r (7) 

where C* = {c^, . . . , c*} is a generating set for £* . A symbol 
node represents a symbol while a check node represents a 
check equation in (|7]l. The check nodes and symbol nodes are 
connected by edges. In fact, we put an edge between j-th 
symbol node and i-th check node if and only if c*^ 7^ 0. A 
Tanner graph for a lattice A is then the Tanner graph of its 
label code (based on construction on £*). 

III. LDPC LATTICES 

There exist many ways to construct lattices based on 
codes |5|. Here we mention two of them. The first one is 
Construction A and the other one is Construction D'. 

A. Construction A versus Construction D ' 

Assume that C is a linear code over Zp where p is a prime 
number, i.e. C C Z^. Let dmm denotes the minimum distance 
of C. A lattice A constructed based on Construction A [5J can 
be derived from C by: 

A = pZ" + e (C) , (8) 

where e: — > is the embedding function which sends 
a vector in Z^ to its real version. In this work, we are 
particularly interested in binary codes and lattices with p = 2. 

Construction D' for lattices (|5] is a good tool for lattice 
construction based on LDPC codes. Let Cq I) Ci 3 • • • 3 Cq 



be a set of nested linear block codes, where Ce [n, fc^, df^^jj^] , 
for < £ < a. Let {hi, . . . , h„} be a basis for F2, where the 
code Ce is formed by the re = n — ke parity check vectors 
hi, ... , hrg. If we consider vectors hj, for I < j < n, as real 
vectors with elements or 1 in M". The new lattice A includes 
all vectors x e Z" satisfying the following modular equations 
hj • X = (mod 2^+1), for < ^ < a and Va^e-i + 1 < .? < 
ra-e- The number a + 1 is called the level of the construction. 

By multiplying the above congruences by proper powers of 
2, one restates Construction D' lfT2l . Indeed, x is in A if and 
only if Hx^ = (mod 2'^+i) where 



H= [hi, 



,2'^h. 



,2'^h. 



(9) 



Then H is called the parity check matrix of the lattice A. The 
lattice A built by Construction D' and associated to this H is 
called an {a + l)-level lattice. Hence the Tanner graph of these 
lattices can be constructed based on their parity check matrices 
H. It can be shown |12| that the volume of an (a + 1) -level 
lattice A is 

det(A) = 22^°=«''^ (10) 
Also the minimum distance of A satisfies the following bounds 



(11) 



A lattice is regular if its Tanner graph, or its corresponding 
parity check matrix H, be (d^, dc) -regular. In other words, 
the number of nonzero elements in all rows and columns are 
the same and are equal to dg and dc, respectively. Now the 
following example reveals the nature of a lattice constructed 
based on Construction D'. 

Example 1: Let a = 2 and Cq, Ci and C2 are 
three nested codes. The sets {0110}, {QUO, 1010} and 
{0110, 1010, 1100}, are the basis for their dual codes respec- 
tively. Therefore, we have hi = (0,1,1,0), ha ^ (1,0,1,0) 
and ha = (1,1,0,0). Let 



H 
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Then x e Z" is in A if and only if Hx^ = (mod 2^). This 
is a 3-level, (2, 2)-regular lattice. 

A lattice A formed by Construction D' is called regular LDP C 
lattice if its parity check matrix H is a sparse (ds, (ic)-regular 
matrix lfT2l . It is easy to see that if the nested codes Ci are 
LDPC codes then the obtained lattice is a LDPC lattice. 

The Extended Progressive Edge Growth Algorithm (E- 
PEG) was introduced and employed to build regular bipartite 
graph 1 12] (regular LDPC lattices). For any < ^ < a, 
denotes the degree of the symbol node s in the parity check 
matrix of Ci. Also for every check node c, the notation dc 
refers to the degree of that node. Suppose that for Q < I < a, 
dl = i + 2 for every symbol node s and select a divisor 
dc — 1°-^^ of n such that dc > dg. Also let be such that 

£ + 2 

re — „ , 1 n (12) 



An (a + 1) -level lattice with these parameters is called an 
(a + 2, 2°+i; a + 1) regular LDPC lattice fHl. These classes 
of lattices are denoted by Af^^^. For sufficiently large n, A°+^ 
is an LDPC lattice with coding gain lfT2l 



7(A 



a+l\ 



4°^ 



(a + l)(a + 2) 

— 



Now, the E-PEG algorithm is employed to construct an 
(a + 2, 2"+i; a + 1) regular LDPC lattice. Indeed, the E-PEG 
algorithm gets a, n and re, < £ < a, as its inputs and returns 
an (a+ l)-level Tanner graph. This Tanner graph corresponds 
to an (a + 2, 2°+^; a + 1) regular LDPC lattice. As we see, if 
the initial inputs of the E-PEG are chosen appropriately, then 
a set of regular Tanner graphs and consequently a class of 
regular LDPC lattices would be constructed. An example for 
a (3, 4; 2) regular lattice is given next. 

Example 2: Consider the (3, 4; 2) regular lattice 

"1000100010010000" 
1000010001001000 
0100100001000100 
0100001010000010 
0010001000100100 
0010000100010010 
0001010000100001 
0001000100001001 
2000002200000200 
0200020000220000 
0020200000002020 
0002000022000002 

Here a = 1, dc, = 4 for 1 < i < 12, = 2 and dl = 3. 
Now we disclose the relationship between a 1 -level lattice 
constructed based on Construction D' and a Construction A 
lattice. 

Proposition 1: Let Aq be generated with Construction D' 
with linear code Co and a = 0, then it is equal to a lattice 
Ai constructed following Construction A using the same 
underlying code Cq. 

Proof: Suppose that Aq be the lattice generated by Cq 
based on Construction D' and a = 0. Also let us assume 
that Ai = 2Z" + e (Co) is the lattice constructed following 
Construction A. The vector x G Z" is in Aq if it satisfies 
the congruences hj • x = (mod 2) for 1 < j < tq. Every 
X g Ai can be represented as x = y + c for y S 2Z" and 
c G Cq. Therefore for 1 < :/ < ro we get 

hj • X = hj • (y + c) = (mod 2), 

since y G 2Z" and h^ • c = (mod 2). It concludes that 
Ai C Aq. If X G Aq, then h^ • x = (mod 2) for 1 < j < tq. 
By subtracting an appropriate z G 2Z" from x, we get x— z G 
{0, 1}" or (x - z) G Z^\ Since h^ • (x- z) ee (mod 2) 
for 1 < J < ro, the vector e^^ (x — z) should be in Cq. Thus 
e^^ (x — z) = Co for some Cq G Cq and hence Aq C Ai. It 
turns out that Aq = Ai. ■ 
Definition 1: A 1 -level LDPC lattice A is a lattice con- 
structed based on Construction D' along with binary linear 



LDPC code C as its underlying code. Equivalently, x G Z" is 
in A if Hx^ (mod 2) where H = [hi, ... , h^J^ is the 
parity check matrix of C. 

Based on Proposition [T] it is clear that a 1 -level LDPC lattice 
A can be constructed by Construction A with underlying code 
C and it also can be represented as A ~ 2Z" + e (C). 

B. LDPC lattices are the same as LDLC 

Low-Density Lattice Codes (LDLC) were introduced 
in iflTl . These lattices have attracted attentions recently ||6|. 
An n dimensional LDLC is an n-dimensional lattice with a 
nonsingular generator matrix G for which the parity check 
matrix H = (G^^) is sparse. In this subsection, we show 
that LDPC lattices are LDLC. 

Theorem 1: Let A = A^+i be a regulai- LDPC lattice built 
by Construction D'. Then A is an LDLC. 

Proof: It is shown in 1121 that if we choose Cq such that 
> 2, then Aw, = 2"+iZ and Pw, (A) = Z, for every 
l<i<n. Hence A* = 'E'^ + ^C*. It means that g, = 2°+i 
and in the view of representation (|6j, we have P^^(A) = I„ 
and C"i(A) = i^\n- In the case for which g.i = 2°+^ the 
row vectors of H as in (l9|l constitute the generator for £*. In 
other words, H is the generator matrix of £*, the dual of label 
code £ of A. Since we have constructed LDPC lattices such 
that H is sparse. So a generator matrix for A* can be obtained 
using the method described in III-BI Thus the row vectors of 
H along with row vectors of I„ form a generating set for A*. 
This implies that regular LDPC lattices AJ'+i are LDLC. ■ 

IV. Numerical Analysis of 1 -level LDPC lattices 



A. Generalized min-sum algorithm 

This algorithm is the generalization of min-sum algorithm. 
It can be employed to solve CVP An explanation of this 
algorithm is given in |12|. In order to solve CVP, given an 
input vector y, this algorithm finds the closest lattice vector 
X to y by finding the label codeword of the coset which x 
belongs to. The generalized min-sum is a message passing 
iterative algorithm. In each iteration two steps are performed: 
(1) symbol node operation and (2) check node operation. Also 
in each iteration algorithm does a hard decision. For further 
details, we refer the reader to llT2l . 

B. Simulation results and comments for 1-level LDPC lattices 

It is known that maximum likelihood decoder on the 
AWGN for codes is minimum distance decoder and maximum 
likelihood decoder on the AWGN without restriction for 
lattices is Euclidean minimum distance decoder. Now assume 
that we are going to use the above generalized min-sum 
algorithm in order to decode a 1-level LDPC lattice. In this 
case for 1 < i < n, we have gi — 2, det(Awi) = 2 
and det(P>V;(A)) = L It has two consequences. First, the 
decoding complexity, which is provided in 1(121 . is much 
lower than the case for 2-level ones. Second, with the above 
parameters the initial weights of the generalized min-sum 



algorithm are equal to the squared distance between the i- 
th symbol of the received word and different cosets of the 
i-th coordinate, which are only two cosets. In fact, these 
weights play the same role as the log-UkeUhood costs do in the 
original min-sum algorithm. It means that in the case of 1-level 
LDPC lattice our generalized min-sum algorithm turns to be 
the original min-sum algorithm which uses Euclidean distance 
instead of the negative log-likelihood costs. On the other hand, 
the squared distances between the elements of the label group 
and each coordinate of the received word are bounded above 
by 1, instead of 4 in the 2-level case. Furthermore, preliminary 
results show that by using 1-level LDPC lattices, our decoder 
stops after 10 to 15 iterations. Since we are adding the costs 
during the decoding process, this causes generation of large 
numbers in the case of 2-level lattices in contrast with 1-level 
lattices. These two advantages allow the algorithm runs for 
much higher dimensions. 

The simulation results of 1-level regular LDPC lattices are 
given in Fig. 1. The graphs show that a gain of 1 dB can be 




Fig. 1. Perfomiance of 2000, 6000 and 10000-dimensioiial 1-level regulai' 
LDPC lattices. 



obtained if we employ 1-level LDPC lattices rather than 2-level 
LDPC lattices Iil2| . Also 1-level regular LDPC lattices of sizes 
n = 6000 and n = 10000, at NEP of 10"^ achieves = 
1.15 dB and — 1.1 dB away from capacity respectively. 
At the moment, the performance of these lattices are being 
evaluated under sum-product algorithm as well. The primary 
results show that at n = 10000 these lattices may work as 
close as 0.67 dB at SER 10^^. This is a normal achievement 
in comparison with previous works where LDLC lattices [17j 
of sizes n ~ 10000 and turbo lattices 1161 of size n = 10131, 
at symbol error rate (SER) of 10~^, achieves — 0.8 dB 
and q;2 = 0.5 dB away from capacity respectively. 

V. Concluding remarks 

The 1-level regular LDPC lattices have been analyzed. 
It has been shown that these lattices are equivalent to the 
lattices constructed based on Construction A. It has been also 
established that these lattices are the same as LDLC lattices. 



Experimental results show that they outperform well-known 
LDLC lattices. That is due to the excellent performance of 
employed underlying LDPC codes. In summary, the 1-level 
regular LDPC lattices perform very well among the regular 
LDPC lattices. 

Further researches have to be done in order to find proper 
shaping methods. In that case one can extract good lattice 
codes from LDPC lattices which are appropriate for both 
Rayleigh fading and AWGN channels |3J. 
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